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Abstract: A simple model for QCD is presented, which is able to reproduce the 
meson spectrum at low energy. The model is a Lipkin type model for quarks cou- 
pled to gluons. The basic building blocks are pairs of quark-antiquarks coupled 
to a definite flavor and spin. These pairs are coupled to pairs of gluons with spin 
zero. The multiplicity problem, which dictates that a given experimental state can 
be described in various manners, is removed when a particle-mixing interaction is 
turned on. In this first paper of a series we concentrates on the discussion of me- 
son states at low energy, the so-called zero temperature limit of the theory. The 
treatment of baryonic states is indicated, also. 
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1 Introduction 



QCD is the favored theory of the strong interactions. At low energy, however, 
the description of the hadronic spectrum based on QCD becomes difficult due 
to the non linear structure of the theory. This non-perturbative regimen, in 
contrast to lattice gauge calculations, can be explored by means of schematic 
models. The use of such models is common to other fields of physics where 
the many-body structure of the theory may be explored by introducing ef- 
fective degrees of freedom and their couplings. Nuclear structure physics 
is one of these examples and it is as complicated and involved as the low- 
energy domain of QCD. Like in the case of nuclear structure physics, QCD 
descriptions based on simple models may help in the understanding of basic 
concepts and procedures. The Lipkin model PQ is one of the most famous 
schematic models, and it helped substantially to appreciate the importance 
of pairing two-body interactions as well as the importance of collectivity in 
building the low-energy part of the nuclear spectrum. An extended version 
of the Lipkin model was applied to the description of pion condensates in 
nuclei [2] . A variety of many body techniques have been tested with Lipkin- 
type models PJHj. In [3] some realistic, less schematic, nuclear interactions, 
suitable to describe various nuclear properties, were investigated in this way. 
In Ref. a Lipkin model was applied to describe a system of many quarks. 
As seen in these examples, the predictions of schematic models can be also 
rather rich in their complexities. This fact was shown, for a simple model of 
many gluon systems, in Ref. jE]. 

Until now, the only formalism which can handle QCD from first principles 
is the lattice gauge theory [7j. Particularly, in many gluon systems, a good 
description is obtained without considering finite volume effects jS]. The 
problem with lattice gauge calculations, to treat QCD at low energies, is 
that only the lowest state, and in some cases also the next to the lowest 
state, for a given spin, charge conjugation and parity, can be calculated. 
Lattice theory is numerically quite involved, and the inclusion of quarks 
and antiquarks brings in additional problems which are difficult to solve. 
Effective models of the hadrons, like the MIT model jU], can help to shed 
some light into the structure of QCD at low energy. In the spectrum 
of gluons, as obtained in jHj, was reproduced and the sequence of levels 
explained by simple assumptions. In other works [Tm iTT] many body methods 
were used to describe the spectrum of QCD at low energy. After these 
considerations, it is obvious that it would be nice to have a model which: i) 
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must be able to describe the basic structure of QCD at low and high energy, 
and ii) must be solbable exactly. Probably such a model does not exist, due 
to the complicated structure of QCD. Nevertheless, one can try to construct 
a model which comes as near as possible to QCD. 

The purpose of the present work is to present a model which fulfills the 
above requirements. In Ref. ^2] the most simple version of such a model 
was presented. Like the models mentioned at the beginning, it is based on a 
Lipkin type model and it consists of two levels for the description of quarks 
and antiquarks. These quarks are coupled to a boson level which describes 
gluon pairs coupled to spin zero. The other gluon states are treated as 
spectators. The basic ingredients of the model are the quark-antiquark pairs 
coupled to flavor singlet and spin zero and gluon pairs with spin zero. In 
Ref. ^2] it is shown that the model is able to describe the appearance 
of a quantum phase transition at zero temperature, when the interaction 
is turned on, and a phase transition to the non interacting case at non- 
vanishing temperature. In Ref. ^2] the basic features of the model were 
discussed. Only flavor singlet and spin zero mesons were taken into account. 
The appearance of a Goldstone boson was obtained for a sufficiently strong 
interaction. This state consisted of a meson with negative parity. Also, it 
possesses a very collective nature, i.e. it is a superposition of many particles 
(quarks, antiquarks and gluons) states. The behavior of the model at high 
temperature was discussed, together with some consequences for the Quark- 
Gluon Plasma (QGP) in ^{TQ. 

In the first part of the paper we will introduce the general form of the 
model for the description of meson states. The discussion will concentrate 
on the behavior of the model at low energy, corresponding to the zero tem- 
perature regime of QCD. The study of the high energy behavior and the 
transition to the Quark Gluon Plasma f3J ES wn l be presented in the forth- 
coming paper of the series [T3]. In section 2 the basic ingredients of the 
model are introduced, with the proposition of a Hamilton operator. Because 
of the difficulties to treat fermion pairs exactly, we shall diagonalize them in 
a boson mapping scheme. The basis used to deal with the bosonic images 
of the fermion pairs, and the corresponding matrix elements of the proposed 
Hamiltonian, are given in the same section 2. There we show how to assign 
charge conjugation and G-parity symmetries to the states belonging to the 
basis. In section 3 the model is applied to the description of the low energy 
meson spectrum. Conclusions are drawn in section 4. 
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2 The Model 



As indicated in Ref. the fermion sector is described by a Lipkin type 
model pQ, consisting of two levels, one at energy — Uf and the other at energy 
+Uf (see Fig. 1). This is the Dirac picture for fermions, where antiquarks 
are regarded as holes in the lower level. Alternatively, one also describes 
quarks and antiquarks in the level at +Uf. The quarks and antiquarks are 
coupled to a 1.6 GeV level, which is occupied by gluon pairs with spin zero 
jH], a value which was obtained in the description of a many gluon system 
and which was adjusted to lattice gauge calculations [Hj. In consequence, 
we shall take the energy of the gluon pair state as an externally fixed data. 
The value ujf is fixed to one third of the nucleon mass (0.33GeV). There are 
further gluon states which do not interact with the quarks and antiquarks. 
These states will be treated as spectators and should be taken into account 
in the final spectrum. 

The degeneracy of each fermion level is 2Q, where Q refers to color, flavor 
and eventually other degrees of freedom. If only spin (n s ), flavor (n/) and 
color (n c ) degrees of freedom are considered, the value of Q is given by the 
product 2Q = n c rifn s . For flavor (0,0) and spin pairs, only, the model has 
similarities to the one of Ref. [2] , but with a different interaction. In Ref. [2] 
the pion condensate in nuclei was the dominant phenomena. In the present 
case, nucleons are replaced by quarks and the pions by gluons. The model 
has some similarities to Ref. j3], which is also a Lipkin-type model. There, 
only quarks were considered and the interaction conserves their number. 

For zero temperature and no interactions the lower level is filled by 
fermions. The creation (annihilation) operators of these fermions are ^ 
( c <*(i,o)/o-i^ j n CQ _ anc j con tra- variant notation for the indices. The symbol 
(1, 0)/ refers to the flavor part, where (1, 0) is the SU(3)-flavor notation and 
/ is a short hand notation for the hypercharge Y, the isospin T and its third 
component T z . The index a represents the two spin components ±5, the in- 
dex i — 1 or 2, stands for the upper or lower level and the index a represents 
all remaining degrees of freedom, which are at least 3 because of the color 
degree of freedom (when only color is taken into account, instead of a we 
will use the index c) . Lowering and raising the indices of the operators intro- 
duces a phase, which depends on the convention used [TB], and a change of 
the indices to their conjugate values, i.e., the quantum numbers (1, 0)YTT z a 
change to (0, 1) — YT — T z — a. 

The operators, defined above, contain the relevant degrees of freedom of 
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Figure 1: Schematic representation of the model space. The fermion levels 
are indicated by their energies ±Uf. The gluon-pairs are represented by the 
level at the energy cuj,. 
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QCD, i.e. color, spin and flavor. These basic degrees of freedom appear at 
all energies, no matter how the resulting particles are defined, i.e., either in 
the perturbative or in the non-perturbative regime. In the non-perturbative 
regime one usually denotes them as constituent or effective particles. This is 
mainly due to the difference in the spatial properties, while color, spin and 
flavor have the same meaning as in QCD. Here, the quarks and antiquarks 
are constituent particles at low energy and have little in common (except for 
the quantum numbers mentioned) with the ones at high energy. We shall 
show that a model which contains these basic degrees of freedom and which 
takes into account the dynamic coupling with gluons can describe the main 
characteristics of QCD at low energy. 

The quark and antiquark creation and annihilation operators are given 
in terms of the operators c and 



a afa- — C a/crl) d a fa- — C afa2 

a afa = c af*l^ d ^afa = c afa2 ^ ^ 

which corresponds to the Dirac picture of particles and antiparticles: quarks 
are described by fermions in the upper level and antiquarks by holes in the 
lower level. 

The gluon sector of the model space is described by bosons which rep- 
resent pairs of gluons coupled to spin zero. The energy of a boson state is 
fixed at the value u>b = 1.6 GeV jSj, as mentioned before. 

The quark-antiquark pairs of the model are given by 



W20-22 _ Dt/202 _ „t ^a/20-22 _ _t jfo/20-2 

a a 
W20-2I 13/20-2 _ t n af 2 <T2l — A n af2&2 

a a 

(~tf20- 2 l _ V^-t ^0/20-2! _ „t „o/20"2 

Wl<7il ~ Z^ C a/i^il C - 2^ a af 1 a 1 a 

a a 

f-if'20- 2 2 t „a/ 2 o-22 _ J r f]af2<T2 (n\ 



The first two equations describe the creation and annihilation of quark- 
antiquark pairs. The pairs can be coupled to definite flavor (A, A) = (0, 0) 
or (1,1) and spin S — or 1. We shall write, in this coupling scheme, 
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as the configuration where the lower state is completely filled and the upper 
one is empty. Note, that the vacuum state is not uniquely defined [T2"] . All 
states, which contain only quarks in the upper level and where the lower 
level is completely filled (so that antiquarks are not activated), regardless 
of color, as for example the three quark baryon states, are annihilated by 
B(x,\)f,SM- This property derives from the fact that the operators B(a,a)/,Sm 
contain an antiquark annihilation operator which anticommutes with all the 
quark creation operators. Therefore, the Hilbert space of the model may be 
divided into sectors, each one with a different vacuum state having a given 
baryon number. The one with the baryon number zero is the real particle 
vacuum. 

2.1 Group Theory of the Fermion Part 

From now on, we restrict to 2Q = n s n c nf = 18 with n s = 2, n c = 3 and 
rif = 3, for the spin, color and flavor degrees of freedom, respectively. The 
largest group, whose generators are cl l j- iaii c C2 ^ 2 ' 72 ^ (cj = 1, 2, 3, fa = 1, 2, 3, 
ai — 1,2 and i,j = 1,2), is the £/(4f2) group. One possible group chain for 
the classification of the states, which include the flavor (SUf(3)) and the spin 
(SUs(2)) groups, is given by 



where the irreducible representation (irrep) of U(AQ) is the completely anti- 
symmetric one and N is the number of particles involved. The upper index 
in [h T ] refers to the transposed Young diagram of [h] , where the columns and 
rows are interchanged [TJj- Due to the antisymmetric irrep [1^] of U(AQ) 
the irreps of U(Q/3) and £7(12) are complementary and the irrep of U(Q/3), 
which is for Q = 9 the color group, has maximally three rows 1 7j. In the 
group chain Q no multiplicity labels are indicated. There is a multiplicity 
Pf for (A/,/i/) and ps for the spin S. The color labels (XcPc) are related 
to the hi via Xc = hi — h 2 and pc = h 2 — h 3 . The complete state is given by 



[1*] 
U(4to) 



[h] = [hifafa] [h T ] 
D 17(g) ® £7(12) 

U U 



(\ c ,Pc) SU C (3) (Xf,p f ) SU f (3)®SU s (2) S,M 



(3) 
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N, (A c , fic),Pf(X f , fMf)YTT z , Ps SM > 



(4) 



where Y is the hypercharge, T is the isospin and T z its third component. For 
meson-like states, the color quantum numbers to be considered are (Ac, pc) = 
(0, 0). These states will be located in an elementary volume of about ^fm 3 , 
corresponding to a sphere of radius 1 fm. 

To obtain the values of hi one has to consider all possible partitions 
of TV = hi + h 2 + h 3 , which fixes the color. For colorless states we have 
hi = h 2 = h 3 = h. Each partition of iV appears only once. The irrep [hhh] 
of If(§) = 17(3) (Q = 9) fixes the irrep of 17(12), as indicated in ©. For 
the reduction of the irrep of £7(12) we have written a computer code [TH] . 
which is available to the interested reader. As an example, let us consider 
the If (12) irrep [3 6 6 ] and the two 1/(4) irreps [9 2 2 ] and [9720], where the 
first one contains the state where the lower level is completely filled and the 
upper one empty, and the second irrep is the next highest one which contains 
flavor (0,0) . The first is accompanied by flavor (0, 0)i and the second one 
by (0, 0)i and (1, l)i, where the subindex denotes the multiplicity. The spin 
content of [9 2 2 ] is given by 55 , 1 45 , 2 36 , 3 28 , 4 2 i, 5i 5 , 6i , 7 6 ,8 3 and 9 t . The 
spin content of [9720] is Ogi, lm, 2i 8 g, 3i35 , 4 90 , 554, 627 and 7g. The lowest 
dimensional irrep is [5 2 4 2 ] with the spin content 0i and li. 

2.2 The Boson Mapping 

The explicit construction of the basis states, Eq. (J3J) and the calculation of 
the matrix elements, become very involved, which is in conflict with the idea 
to develop a simple model. A way out of it is to use a boson mapping of the 
pair operators B^ and B and work in the boson model space. 

The basic ingredients of the model are the pair operators, given in (0). 
They can be mapped onto boson operators (T^j 



Ej/20-2 L/20-2 /r\ 

where the operators on the right satisfy the normal boson commutation re- 
lations. 
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- ^/ 3 /2<Wi<W 2 <Wi • (6) 

The exact boson mapping is quite involved, but it can be obtained in general 
[T9[ 1201 EI]- For the sake of this work, it is worth to show that the mapping 
can be performed, indeed. We shall: i) work from the beginning in the boson 
space, ii) define a Hamiltonian which corresponds to the Hamiltonian acting 
in the fermion space, and iii) select a basis. The advantage of working in the 
boson space is the simplification in getting the matrix elements (see below). 
The price to pay is related to the appearance of non-physical states [TH] , as 
we shall discuss later on. 

In order to choose a basis in the boson space, we profit from the fact 
that the basic degrees of freedom are given by the boson creation operators 
b(\,X)fSM = b\fSMi w ith A = or 1 and S = or 1. This gives four possible 
combinations of [A, S]: [0,0], [0,1], [1,0] and [1,1]. Consequently, the total 
Hilbert space is the direct product of a one, three, eight and 24 dimensional 
harmonic oscillators [22]. 

For each harmonic oscillator we can define a seniority basis 



/20-2 tJUo-i 
fieri' /30-3 



^N xs u xs (b[ s ■ b\ s ) 2 Ks«as> , (7) 

where N xs is the number of bosons of type [A, S], v\ S the corresponding se- 
niority and Nn xs v xs is a normalization constant. The seniority is the number 
of ^As-bosons not coupled into pairs. The ci\s contain all other quantum 
numbers for a particular harmonic oscillator. The dot in the factor refers to 
the scalar product. 

The choice of a seniority basis is particularly useful for the calculation 
of the matrix elements of the interaction, which contains expressions of the 
form (b\ s -b xs ), (b\s-b\s) and (b\ s -b\ S ), where the latter is just the number 
operator of the bosons of the type [A, S}. The exact structure of I^as^as > 
is not needed, except for the knowledge of the quantum numbers a\s- 

For the one dimensional harmonic oscillator ([0,0]) the seniority can take 

the values or 1. The state is of the form (bl o ) Noo \0 > = (bJo&oo)^ 2 ^ 
{bl ) Uoo \0 >. For the three dimensional harmonic oscillator the seniority is 
equal to the spin S\s- The explicit expression of the state is given in Ref. 
|22j . The eight dimensional oscillator contributes to flavor only and it is 
discussed in Appendix A. The 24 dimensional oscillator can be found in Ref. 
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jHj, where the color part in Ref. jB] has to be interpreted here as the flavor 
part. In Ref. jB] only singlet states are listed, but the procedure to obtain 
non-singlet flavor states is outlined. 

The parity of each state is given by P = (— 1)^, where N = J2\,s N\s is 
the total number of bosons. Each boson stems from a particle-antiparticle 
pair, which carries negative parity. 

In order to obtain the property under charge conjugation, one has to apply 
the charge conjugation operator C to the pair creation operator B^ X \)f sm- 
The result is (see Appendix B) 



CB Xf SM C 1 - (-1) B- x j SM , (8) 

where A = (A, —fi), f = —Y, T, —T z . From this it is clear that only states 
with Y = 0, T z = 0, and p = 0, can have a definite C-parity. In Eq. (jHJ) we 
make use of the application of the operator C which interchanges quark and 
antiquark operators (a* <-> S) and inverts the magnetic quantum numbers 
(Yi — > —Y{ and Ti z — > —Ti Z ) of flavor and of color only 1 . 
For products of two pair creation operator we obtain 



C 



f,M 



c- 1 



\Sl + S2— A — H+Pr, 



B 



B 



p(/i,A) 

f,M 



(9) 



where p is the multiplicity label of (A,/i) in the product (A 1; Ai) ® (A 2 , A 2 ). 
The symbol p ma x denotes the maximal value of p. The phase convention of 
Ref. was used. The symbol ® denotes the combined product in SUf(3) 
and Si/a (2). 

In analogy, the action of the charge conjugation on a product of three 
pair operators can be obtained: 



C 



B 



B 



^2,S2 



P12(Al2,/il2),S'l2 „l 

® B J 



^ ^Sl+S2+53 — A — fJ,+pi2,max-pi2+Pm.ax— P ~Q 



Al,5l 



B 



t 

M,S2 



ip(a,a0,s 

C 1 

f,M 

Pl2(A t 12,Ai2),5l2 



B l,s 3 



-\p(p,,\ 

f,M 



1 After here, for the sake of notation, we shall indicate charge conjugate states with a 
bar on the index f 
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where pi 2 is the multiplicity of (Ai 2 , /ii 2 ) in the product of (Ai, Ai)®(A 2 , A 2 ), p 
is the multiplicity of the total irrep in the last coupling of the above equation, 
and pi2,max is the maximal value of pi 2 . 

The procedure outlines here can be used in a recursive way for more 
involved coupling schemes. For our purpose it is sufficient to go up to three 
pairs, which will be the dominant structure at low energy. 

For the G-parity the additional rotation in the isospin space has to be ap- 
plied, which changes Tj )2 to — T^ z |24j . For a polynomial in the pair operators 
this gives an additional phase (— 1) T , where T is the total isospin [21]. 

The same phase properties under C- and G-parity transformation have 
to be valid for the mapped boson operators b' X j SM . 

In a seniority basis, the matrix elements are easily obtained, and they are 
written 



< N xs + 2v xs a x s\(bl s ■ b\ s )\N xs v xs axs > 

< N xs ~ 2vxsaxs\(b xs • b xs )\Nxs^xsaxs > 

< Nxsv\s<x\s\(b\s ■ bxs)\NxsV\saixs > 



= \J (N xs - v\s + 2) (N xs + v\s + d xs ) 
= \l (Nxs - v\s) (Nxs + v\s + dxs ~ 2) 
= Nxs , 

(11) 



where dxs is 1> 3, 8 or 24 for the case of the one, three, eight or 24 dimensional 
harmonic oscillator. As a short hand notation we will use instead of (b\ s -b\ s ) 
the expression (^ 5 ) 2 , and similarly for the other products, (bxs) 2 and b xs bxs- 
As a Hamiltonian, invariant under rotation, charge conjugation and G- 
parity, we propose 



// 



2u f n f + uj b n h + 



XS 



[b{ s f + 2b[ s b xs + (bxsf 



^)b + 

2ir 



& f a-g) [(b{ s )" + 2b{ s bxs + (bxs) 



(12) 



Due to symmetry arguments, the interaction strength V x $ is the same for the 
two last lines in Eq. ([12)1 . The term (b\ s ) 2 ((bxs) 2 ) describes the creation 
(annihilation) of two quark-antiquark pairs with the simultaneous creation 
or annihilation of a gluon pair. The term b xs bxs, in Eq. (|12|l . describes 
the scattering of a fermion pair with the emission or annihilation of a gluon 
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pair. All processes can be depicted by a Feynman graph and all graphs 
can be obtained from any other one by an appropriate interchange of lines. 
Because the strength V\s should be, basically, invariant under the exchange 
of lines, we shall use the same interaction strength for all channels, as a first 
approximation. The terms which appear in Eq. (|T2*|) originate in the normal 
product of : (b^ s + bf S ) 2 :, where the square implies a scalar product. The 

factor (1 — 2n ) represents a cutoff which can be traced back to the boson 
mapping of the fermion pairs with flavor (0,0) and spin 0. This term simulates 
the effect of an exact boson mapping I2H1 E] and it is responsible for the 
disappearance of the interaction when the number of pairs reaches 2Q. In 
other words, this cutoff term simulates the Pauli principle which does not 
allow more than 2Q pairs. 

The Hamiltonian (|12j) is the most simple form we can think of and it 
contains only four parameters (the values of V\s). The value of uif is fixed 
to 0.33 GeV, which is about | of the mass of a nucleon. The most notorious 
difficulty, associated to the use of the boson mapping, lies in the Hilbert 
space of the boson operators. It is larger than the Hilbert space of the 
fermion pairs. In some situations one can identify the source of the spurious 
dimensions. If, for example, only flavor (0,0) and spin pairs are taken into 
account, the relevant group structure is £7(40) D £7(20) ®U{2). The irrep of 
£7(40) has to be antisymmetric, which implies that the irreps of £7(20) and 
U(2) have to be complementary, i.e. if U(2) is given by a Young diagram 
with two rows, the one of £7(20) has to be the adjoint, which is obtained 
by interchanging rows and columns ^7j. The upper limit, up to which no 
spurious states appear, is 20 because £7(20) allows 20 rows in the Young 
diagram. This is also the maximum number of pairs allowed, i.e. for this 
case no un-physical states occur. If flavor values (0,0) and (1,1) and spin 
are used, only, we have £7 (40) D £7(^p) <8> U (3) and up to pairs there is no 
problem with respect to the appearance of un-physical states. This implies 
that states with explicit flavor will present un-physical states only for large 
number of bosons. If flavor (0,0) and spin and 1 are considered, we have 
£7 (40) D £7 (0)®£7(4), which has as an upper limit the number O up to which 
no spurious states appear. Finally, for all pairs, i.e. flavor (0,0), (1,1) and 
spin 0, 1 the relevant group chain is £7(40) D £7(§) <g> £7(12) and the upper 
limit is j. This gives us a hint about the group sequence where un-physical 
states may appear. The upper limit up to which all states are physical is 
lowered in the sequence where bosons with more degrees of freedom appear. 
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There are states which can be described either by one or the other type of 
bosons or even by a combination of bosons. For example, when both levels, 
the lower and the upper ones, are filled there is only one allowed state, which 
is the flavor singlet with spin zero. However, all types of boson pairs can 
describe it, e.g. when the number of bosons coupled to flavor singlet and 
spin zero is equal to 2Q. In view of these considerations, and for the one, 
three, eight and twenty-four dimensional harmonic oscillator basis, we have 
introduced the limits 2f2, fi, and p, respectively. The higher non physical 
states do not play an essential role because, as shown below, the dominant 
contribution at low energy comes from configurations with a small number of 
quark- ant iquark pairs ^3] • By working with these dimensional cut-off values 
the influence of non-physical states is minimized. Also, for each case, the 
total number of bosons is restricted to < 2Q. For a reasonable interaction 
strength, however, the dominant contribution comes from a small number of 
bosons. In such cases, the number of un-physical states is small and they do 
not influence much the result. The dimensional cut-off in the Hilbert space 
has the advantage that most un-physical states are excluded. In principle, 
one can eliminate the spurious states by applying another, more involved, 
procedure. For that one has to reduce the irrep of U(12) to the flavor and 
spin groups, as done in the last section. This gives us the allowed content 
of flavor and spin for a given irrep of £7(12). The matching condition, i.e. 
by comparing for a given number of pairs the spin and flavor content on the 
boson side to the one on the fermion side, eliminates un-physical states. If 
on the boson and fermion sides, for a given flavor and spin, the number of 
states are equal, all states in the model space are taken into account. This is 
the case of low lying basis states. A simple counting procedure can be used 
for other situations. If there are, for a given flavor and spin, more states in 
the boson space than in the fermion space, one can not decide easily which 
combination is allowed. However, one can reduce the number of states of 
the model space to the same dimension as the one of the fermion space. As 
a rule one can first eliminate the states which contain most of the bosons 
with the largest degree of freedom, i.e. which are of the type [1,1], and in 
this way proceed, if necessary, until only states with flavor (0,0) and spin 
zero bosons are left. At least, the proposed procedure eliminates most of the 
spurious states. The error made can be absorbed in the parameters of the 
model, a general practice in dealing with phenomenological models, because 
in the end the correct number of degrees of freedom (dimension of the Hilbert 
space) dominates in a successful description of the spectrum. The idea of the 
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proposed procedure is not new and it was used in another context by J. Cseh 
et al. • 

3 Application to the meson spectrum 

The Hamiltonian (j!2|) commutes with the isospin operator and it does not 
depend explicitly on the hypercharge Y. As a consequence, all states which 
belong to the same flavor irrep are degenerate. In principle, we can add terms 
proportional to T 2 , Y and Y 2 in order to lift the degeneracy. These terms 
will add new parameters to the four already present (ub is fixed, as in Ref. 
jHl). Also, a flavor mixing term could be added, as suggested by the 77-77' 
mixing [26 . In order to simplify the discussion, we shall first ignore these 
additional interaction terms. 

In fitting the spectrum of the mesons we will use, as an experimental 
input for each multiplet, only the state with T = and Y = 0. For an octet 
all states have the same energy as the isospin singlet and hypercharge zero 
state. Because later on we shall take into account flavor mixing interactions 
too, the position of the singlet and octet state are not fixed at the measured 
energy values but at the values obtained when the flavor mixing interactions 
are switched off [21] • The mixing angle is introduced for two multiplets: 
the (1,1) 0~ and (0,0) 0~ irrep, containing the pions and the 77, 77', and the 
(1,1) 1~, (0,0) 1~, containing the u, (ft and p particles. The mixing angles 
are, respectively, —23.7° and 35.3°, For other multiplets one assumes 

that the mixing angle is zero, because of missing data, and because of the 
smallness of the energy splitting between members, as compared to the energy 
splitting within the multiplet which contains the pions or the p mesons. The 
uncorrected masses for, e.g., the octet (before flavor mixing) are wig = 615 
MeV in the first case (see notation of Ref. [23) an d 7n 8 = 940 MeV for the 
second case (see also Table 1). 

In Table 1 we show the states used in the fit. Their flavor, spin and parity 
are indicated together with the experimental values. In total, to perform 
the fit, we have considered thirteen states with spin zero and one in the 
four parameter fitting procedure. All other states are predicted, particularly 
those with spin 2 and 3. 

In Figure 2 we give the spectrum for spin and 1 meson states without 
any particle number changing interaction. On each side of a level the flavor 
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Figure 2: The meson spectrum for spin and 1 states, for the case of no 
interaction. The value mf — 0.33 GeV was used. Note the large multiplicity 
appearing at already low energies. 
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77(1760) 


(1,1) 


0" 


1.671 


1.760 


cj(782) 


(0,0) 
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cj(1420) 
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\- 


1.389 


1.420 


cj(1600) 


(1,1) 


\- 
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Table 1: States used in the fit. The particles are listed in the first column, 
and their transformation properties in flavor and spin are shown in the second 
and third columns. Note that, for the particles in the first (0,0), (1,1) 0~ an d 
(0,0), (1,1) 1~ irreps, we are listing the value of the masses without flavor 
mixing (they are marked by an asterisk). The experimental data are taken 
from jZ7] 

quantum numbers and its degeneracy are indicated. This serves to illustrate 
that the multiplicity at energies lower than 2 GeV is already very large. This 
is a consequence of the various manners in which the same set of quantum 
numbers can be obtained, for a given configuration, when many quarks, anti- 
quarks and gluons are considered. This is known as the multiplicity problem. 
The result of the best fit values, obtained after the interaction is turned on, 
is given in Figures 3-6 for spin 0, 1, 2 and 3 respectively. Only states which 
correspond to non exotic parity/charge conjugation quantum numbers are 
shown. Most of them appear above 2 GeV and some can be deduced from 
the gluon spectrum published in Ref. jB]. In Figures 3-6 each theoretical 
spectrum is compared to the experimental one. On the right hand side of 
each level the theoretical interpretation in terms of flavor and the multiplic- 
ity is indicated. On the left hand side of each spin (J PC ) the experimental 
information is given, taken from the particle data group ^7\. The energy of 
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these states, appearing in the summary table of Ref. j2Z], is given in boxes 
and the experimental error is reflected by the size of the box. If the error is 
very small, the box is replaced by a line. States which are not in the sum- 
mary table of [23 are indicated by dashed boxes (lines). Only states which 
correspond to isospin singlet and hypercharge zero, after having corrected 
for the isospin mixing, are listed. 



Figure 3: The meson spectrum for spin as obtained from the fit to exper- 
imental data [27| . 

Note that very few states have a multiplicity. Most states were pushed 
upwards due to the interaction. This is an effect of the interaction, because it 
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Figure 4: The meson spectrum for spin 1 as obtained in a fit to experimental 
data HZ!- 
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Figure 5: The meson spectrum for spin 2, obtained with the paremeters 
fixed by the fitting procedure. Experimental data are from fI7\ . 
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Figure 6: The meson spectrum for spin 3. Experimental data are from 
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changes the number of particles and relates the quark-antiquark sector with 
the gluons. Models with a particle conserving interaction will always present 
the multiplicity problem. Thus, the particle mixing interaction is essential 
to remove the multiplicity problem. 

The spin 2 and 3 states where not adjusted. The theoretical results seems 
to agree with the data, because the states are predicted at the correct energy 
domain. The density of states for a given flavor also seems to be reproduced. 



particle 


Etheo 




r 


< n w > 


<n q > 


<n g > 


vacuum 


0.0 


(0,0) 


0+ 


3.118 


3.177 


1.705 


/ (400-1200) 


0.656 


(0,0) 


0+ 


0.457 


0.471 


0.321 


/o(980) 


0.797 


(1,1) 


0+ 


3.781 


3.832 


1.495 


/i(1420) 


1.445 


(0,0) 


1+ 


2.392 


3.434 


0.902 


/ 2 (1270) 


1.363 


(1,1) 


1+ 


2.464 


3.519 


0.993 


7/(958) 


0.885 


(0,0) 


0- 


2.509 


3.562 


1.292 


7^(1440) 


1.379 


(0,0) 


0- 


0.773 


1.790 


0.444 


77(541) 


0.602 


(1,1) 


0- 


2.711 


2.766 


1.163 


77(1295) 


1.428 


(1,1) 


0- 


1.611 


1.638 


0.531 


t?(1760) 


1.671 


(1,1) 


0- 


3.535 


4.581 


1.254 


u(782) 


0.851 


(0,0) 


\- 


2.563 


3.621 


1.341 


0(1020) 


0.943 


(1,1) 


\- 


2.394 


3.438 


1.198 


w(1420) 


1.389 


(1,1) 


\- 


0.853 


1.870 


0.468 


cj(1600) 


1.639 


(1,1) 


\- 


3.546 


4.597 


1.206 



Table 2: Particle content for selected states. In columns we indicate the 
theoretical energy (Etheo), the flavor ((A/,///)), spin J and parity (ir), expec- 
tation value of the boson pairs in the channel (1,1) 0~, (< 77,10 >), expectation 
value of the total number of quark-antiquark pairs (< n q >) and the total 
number of gluon pairs (< n g >) with spin 0. 

In Table 2 we show the quark-antiquark pair and gluon pair contents of some 
selected states. The total number of quark-antiquark pairs is denoted by 
< n q >, where the symbol < ... > indicates the expectation value of this 
number. The quantity < > gives the average number of boson pairs of 
the type [i,j], while < n g > is the expectation value of the number of gluon 
pairs with spin zero. The total number of gluons is twice < n g >. 
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The structure of the calculated ground state (physical vacuum) is an 
interesting piece of information about the model predictive power. The cal- 
culated value of the ground state energy is equal to -0.726 GeV. The physical 
vacuum state contains about 3.1 quark-antiquark pairs of type [1,0], and the 
other configurations contribute with 0.06 pairs. The dominance of the [1,0] 
quark-antiquark pairs is consistent with the strengths of the interactions. 
The parameters, obtained in the fit, are : m/ = 0.33 GeV, Vo,o = 0.0337 
GeV, Vo,i = 0.0422 GeV, V lj0 = 0.1573 GeV and V ltl = 0.0177 GeV. The 
channel [1,0] is clearly the strongest one. Also, the number of gluons is no- 
ticeable, i.e 1.7 pairs, which correspond to more than three gluons contained 
in the elementary hadron volume. This indicates a collective behavior of 
QCD states at low energy. Indeed, the pion state (corresponding to 77(541) 
in Table 2) contains about 2.7 pairs of the type [1,0] while the rest contributes 
with about 0.06 pairs, as in the ground state. The number of gluon pairs in 
the pion, 1.2 pairs or 2.4 gluons, is similar to the number of gluon pairs in the 
ground state. They constitute about 30 per cent of the particle content. In 
brief, the calculated states contain a large number of quark-antiquark pairs 
and gluons. Roughly speaking, no single state can be described approxi- 
mately by a pure quark-antiquark pair. Note that in this respect theory and 
experiment do agree, in spite of the simplicity of the model. We think that 
the complex structure of the meson spectrum is described qualitatively by 
our model, as well as several other features, like the position of the first states 
with spin 2 and 3, and the density of states with flavor (0,0) or (1,1) for each 
spin and parity, charge conjugation number. 

Concerning baryons, we have to include them yet in the model. Up to 
now they are described as spectators, i.e. without having an explicit coupling 
to the quark, antiquark, and gluon sea. For that, further interaction terms 
should be introduced. For example, one can introduce the interaction 

nfl,(o,i)o(b t + b) + tifl,( 2) o)i(b t + b) , (13) 

where nD,{x,n) s is the number operator of a Di-quark coupled to flavor (A, //) 
and spin S. This is analogous to the above ansatz of the Hamiltonian. The 
product of two pair creation operators of Di-quarks can not appear because 
this would mix the baryon number. The interaction in (|13|) is a direct ex- 
tension from the pair operator interaction of the former Hamiltonian. The 
first term in ()13|) acts only on states like the nucleon octet and the last one 
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on particles like the baryon decouplet. The residual interaction will mix the 
number of gluons, and, mainly, increase the gluon content in the baryons, 
from a gluon content of about 30 percent to, say, 50 percent. The inclusion 
of baryons will be reported in another work [2E] ■ 

4 Conclusions 

In the present paper we have advanced a schematic model of QCD, based on 
a Lipkin-type model for fermions and an interaction to one gluon pair states. 
We have discussed the low energy structure of the model. The Hamiltonian 
is composed by a diagonal, particle conserving part, and an interaction which 
couples the quark-antiquark pairs to the gluons and changes the number of 
particles. The model contains only 4 parameters which were adjusted to 
reproduce 13 observed meson states with spin and 1. After fixing these 
parameters, we have predicted the remaining part of the spectrum. The 
complex structure of the meson spectrum was qualitatively reproduced by 
our results. 

Due to the schematic nature of the model one cannot expect to be able 
to reproduce all details of the low energy meson spectrum. However, the 
results are in qualitative agreement with data, a fact that shows the validity 
of the model as a toy model for QCD. Baryons where not considered but 
the extension to this sector was briefly indicated. The baryons would corre- 
spond to states where three extra quarks are added in the upper level. The 
corresponding operators will commute with the boson pair operators and an 
interaction of baryon states with the quark-antiquark sea should be included 
in the model. 

We have found that the inclusion of particle mixing interactions turns out 
to be essential in order to remove the multiplicity problem encountered in 
other models, when states with many quarks and antiquarks are considered. 
This particle changing interaction also introduces ground state correlations 
resulting in many quarks, antiquarks and gluons configurations in the states. 
It produces a large contribution of the gluons and the total spin is not a 
simple product of a quark-antiquark state but of many quarks, antiquarks 
and gluons. This illustrates the fact that, even at low energy, the structure 
of the hadron states is by no means as simple as suggested by earlier particle 
conservation models. It also shows that phenomenological potentials, which 
simulate the presence of gluons in a pure quark model, cannot resolve the 



23 



problem of multiplicity. 
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Appendix A: The eight dimensional oscillator 

The reduction of the eight dimensional oscillator group, £7(8), to the flavor 
group SUf(3) is discussed in Ref. jSJ. As an intermediary group, between 
£7(8) and ££7/(3) one can use the 50(8) group. Though, in Ref. [S] only the 
reduction to flavor singlet groups is listed, the general procedure is outlined. 
Programs are available on request and the procedure has been published 
elsewhere |3Ti] . 

The generators of the U(8) group are given by bL-^^ oo^(i,i)/2,oo, where 
the zeros refer to zero spin and its projection. Therefore, these bosons can 
only contribute to the flavor content. In Table 3 we list the flavor content of 
up to four bosons of the type bL i)/ o- As one can see ' the multiplicity raises 
especially for the (1,1) flavor irrep. With the help of the 5*0(8) group one 
can further reduce the multiplicity. For our purpose this is not necessary. 

Appendix B: Parity, Charge Conjugation and 
G-Parity 

The charge conjugation operator acts as follows on the quark and antiquark 
creation operators 
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Table 3: The first column gives the Young diagram of the U(8) group, the 
second column the irrep of the flavor group SUf(3) and the third column 
gives the multiplicity of the flavor irrep. Only the flavor irreps (0,0) and 
(1,1) are listed. 
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Cal fa C-' = dl fa 

Cdl^C- 1 = a\ Ja , (14) 

where the transforms in color and flavor as a (1,0) SU(3) irrep, while the 
d^ transform as (0,1). If one of the color or flavor index is raised then the d) 
transform as (1,0). The "c" and "/" refer to the reflection in the magnetic 
quantum numbers of SU{3). I.e. / stands for Y, T and T z and / for —Y, T 
and — T z and similar for " c" . 

With this, the action of the charge conjugation operator on a quark- 
antiquark pair is given by 



CB^C- 1 = CE c/lMff2 4 fiai d\; a2 < (1,0)/ 1 ,(0,1)/ 2 |(A 1 ,A 1 )/ >i 

= E cflf2 a^4 ha a}^ < (l,0)/ 1 ,(0,l)/ 2 |(A 1 ,A 1 )/>i {\a u \a 2 \SM) 
= -Ecf lh ^ 2 a\- 2 c a dl hai < (l,0)/i,(0,l)/ 2 |(A 1 ,A 1 )/ >i {\^\a 2 \SM) 
= -(-l) 2Al+1 - 5 E c/l/2CTlCT2 a} 2 c a dl fiai < (l,0)/ a ,(0,l)/i|(Ai,Ai)/> 1 {\a u \a 2 \SM) 

= (- 1 ) SB h,s 1 M . (15) 

where we made use of the properties of the SU (2) [HI] and SU(3) [23\ Clebsch- 
Gordan coefficients. The subindex 1 in the SU(3) Clebsch-Gordan coefficient 
indicates a multiplicity of one j2H| • 

For the product of two pair operators we have 
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where p is the multiplicity index of (A,//) in the product (A 1; Ai) ® (A 2 , A 2 ). 
For the product of three pairs we have 
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with the use of the notation of Ref. [23] for the SU(3) Clebsch-Gordan 
coefficients and their symmetry properties. 
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